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Abstract. We present a computer simulation model for the Hanbury Brown-Twiss ex- 
periment that is entirely particle-based and reproduces the results of wave theory. The 
model is solely based on experimental facts, satisfies Einstein's criterion of local causal- 
ity and does not require knowledge of the solution of a wave equation. The simulation 
model is fully consistent with earlier work and provides another demonstration that 
it is possible to give a particle-only description of wave phenomena, rendering the 
concept of wave-particle duality superfluous. 
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1 Introduction 

Computer simulation is widely regarded as complementary to theory and experiment [1]. 
Usually, the fundamental theories of physics provide the framework to formulate a math- 
ematical model of the observed phenomenon, often in terms of differential equations. 
Solving these equations analytically is a task that is often prohibitive but usually it is 
possible to study the model by computer simulation. Experience has shown that com- 
puter simulation is a very powerful approach to study a wide variety of physical phe- 
nomena. However, recent advances in nanotechnology are paving the way to prepare, 
manipulate, couple and measure single microscopic systems and the interpretation of the 
results of such experiments requires a theory that allows us to construct processes that 
describe the individual events that are being observed. Such a theory does not yet exist. 
Indeed, although quantum theory (QT) provides a recipe to compute the frequencies for 
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observing events, it does not describe individual events, such as the arrival of a single 
electron at a particular position on the detection screen [2-5]. Thus, we face the situation 
that we cannot rely on an established physical theory to build a simulation model for the 
individual processes that we observe in real experiments. Of course, we could simply 
use pseudo-random numbers to generate events according to the probability distribu- 
tion that is obtained by solving the Schrodinger equation. However, that is not what the 
statement "QT does not describe individual events" means. What it means is that QT 
tells us nothing about the underlying processes that give rise to the frequencies of events 
observed after many of these events have been recorded. Therefore, in order to gain a 
deeper understanding in the processes that cause the observed event-based phenomena, 
it is necessary to model these processes on the level of individual events without using 
QT. The challenge therefore is to find algorithms that simulate, event-by-event, the ex- 
perimental observations that, for instance, interference patterns appear only after a large 
number of individual events have been recorded by the detector [4, 6], without first solv- 
ing the Schrodinger equation. 

In this paper, we leave the conventional line-of-thought, postulating that it is fun- 
damentally impossible to give a logically consistent description of the experimental re- 
sults in terms of causal processes of individual events. In other words, we reject the 
dogma that there is no explanation that goes beyond the quantum theoretical description 
in terms of averages over many events and search for an explanation of the experimen- 
tal facts in terms of elementary, particle-like processes. It is not uncommon to find in 
the recent literature, statements that it is impossible to simulate quantum phenomena by 
classical processes. Such statements are thought to be a direct consequence of Bell's the- 
orem [7] but are in conflict with other work that has pointed out the irrelevance of Bell's 
theorem [8-31]. This conclusion is supported by several explicit examples that prove 
that it is possible to construct algorithms that satisfy Einstein's criterion for locality and 
causality, yet reproduce exactly the two-particle correlations of a quantum system in the 
singlet state, without invoking any concept of QT [32-37]. It is therefore an established 
fact that purely classical processes can produce the correlations that are characteristic for 
a quantum system in an entangled state, proving that from the viewpoint of simulat- 
ing quantum phenomena on a digital computer, Bell's no-go theorem is of no relevance 
whatsoever. 

This present paper builds on earlier work [32-35,37-45] that demonstrates that quan- 
tum phenomena can be simulated on the level of individual events without first solving 
a wave equation and even invoking concepts of QT, wave theory or probability theory. 
Specifically, we have demonstrated that it is possible to construct event-by-event proce- 
ses, that reproduce the results of QT for single-photon beam-splitter and Mach-Zehnder 
interferometer experiments [6], Einstein-Podolsky-Rosen-Bohm experiments with pho- 
tons [46-48], Wheeler's delayed-choice experiment with single photons [49], quantum 
eraser experiments with photons [50], double-slit and two-beam single-photon interfer- 
ence, quantum cryptography protocols, and universal quantum computation [40,41]. Ac- 
cording to the theory of quantum computation, the latter proves that at least in princi- 
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Figure 1: Schematic picture of a HBT experiment [59]. Top: Source. Coherent light, generated by a YAG 
laser, is sent through the Gan prism GPl,a single slit S, a beam splitter (a CaCC>3 crystal), an electro-optic 
modulator (EOM) and another Gan prism GP2 to produce two beams A and B as if they would have emerged 
from a double slit separated by 1.3 mm [59]. The EOM is switched rapidly to destroy the first-order coherence 
between beams A and B. Bottom: The interferometer consists of two beam splitters BS1 and BS2 and phase 
shifters <p^ n and (pg n (« = 1,2,3). Light intensity is measured by the three detectors D\, D2 and D3. 



pie, we can construct particle-like, event-by-event processes that can simulate any quan- 
tum system [51]. Some interactive demonstration programs are available for down- 
load [52-54]. 

In this paper, we extend the range of applications of the event-based simulation ap- 
proach by demonstrating that the event-based algorithms, used in our previous work, 
can be re-used, without modification, to build a particle-only simulation model for an- 
other fundamental physics experiment, the Hanbury Brown-Twiss (HBT) experiment [55]. 
The HBT effect refers to a variety of correlation and anti-correlation effects in the inten- 
sities received by two or more detectors from a beam of particles [56-58]. According 
to common lore, when a HBT experiment is performed using single-particle detectors, 
the HBT effect is attributed to the wave-particle duality of the beam. In this paper, we 
present a particle-only model of the HBT effect, demonstrating that it is possible to con- 
struct causal, particle-like processes that describe the experimental facts without invok- 
ing concepts of QT. 

As a concrete realization, we consider a recent HBT experiment [59], a schematic pic- 
ture of which is shown in Fig. [l] A radiation source, a frequency doubled Q-switched 
Nd:YAG laser with wavelength 532nm, is used. The coherent light from this source is 
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split by a beam splitter. The electro-optical modulator (EOM) erases the first-order inter- 
ference of the light [59]. The two beams that emerge are labeled A and £>, see Fig.QJtop). 
Then, the two beams are sent to the three detectors through two beam splitters (BS), see 
Fig. [Tfbottom). After measuring the coincidences of three detectors by means of a triple 
coincidence circuit (TCC), the third-order intensity interference pattern is observed [59]. 

The purpose of this paper is to demonstrate that one can construct a simulation model 
of this experiment that 

• is a one-to-one copy of the experimental setup such that each device in the real 
experiment has a counterpart in the simulation algorithm 

• is event-based and satisfies elementary physical (Einstein's) requirements of local 
causility 

• reproduces the results of wave theory by means of particles only. 

The structure of the paper is as follows. In Section 2, we briefly review the wave the- 
ory of second and third-order coherence. The simulation model is described in Section 3. 
Section 4 presents our simulation results and a discussion thereof. Our conclusions are 
given in Section 5. 



2 Wave theory 

Conceptually, the experiment of Fig. [T] can be viewed as a double-slit type experiment 
with three detectors, as shown in Fig. [2] Assume that source A emits coherent light with 
amplitude a and that source B emits coherent light with amplitude B. Thus, according to 
the superposition principle, the total amplitude falling on the n-th detector (n = 1,2,3) is 

a n = aj***+pj***, (2.1) 

where (pAn (<pBn) is the accumulated phase of the photon travelling from source A (£>) to 
the n-th detector. While the intensity is 

I n = \a n \ 2 = I A + I B +2Rea J 6V*» / (2.2) 

where Ia = \oc\ 2 , Ib = \B\ 2 , and (p„ = (pAn — ^Bn- If the relative phase of a. and B is fixed, 
Eq. (|2.2[) predicts that interference fringes will be observed. If there is no correlation 
between the phases of a and B, there are no interference fringes because 

(In) = (lA) + (h). (2.3) 
On the other hand, the product of the intensities is given by 

I n I m = \a n a m \ 2 = \a 2 e i ^ A " + ^ + B 2 e i ^ B " + ^ 



Figure 2: Schematic picture of third order intensity correlation. Photons emitted from sources A and B are 
registered by three detectors D\, D2 and D3. (p An and (pg n (w = l,2,3) are the phases accumulated during their 
flight from sources A or B to the w-th detector. 



and after averaging over the uncorrelated phases of a and f5, we find 

g£J = (I„/ m ) = (^^) + (J B I B ) + (^I B )|e^»+^)+ e ^ + ^)| 2 

= (/i) + ai)+2(J A J B )(l+cos^ w ) (2.5) 

where nm = <j> n — <pm and n,m = 1,2,3. According to Eq. (|2.5|) the intensity-intensity corre- 
lation will exhibit interference fringes, a manifestation of the so-called Hanbury Brown- 
Twiss effect. It is convenient to introduce the normalized, dimensionless, correlation by 

r (2) 

(2) "ran ,n, r \ 

where (I n ) = (I m ) = (I A ) + (I B ). Assuming that the sources A and B have the same statistics 
and the same average intensities, we have I a = h arid obtain 

& =g {2) (l + ^cos</> nw ) , (2.7) 

where g^> = ( ) / ( 7^ ) 2 is the second-order normalized intensity autocorrelation func- 
tion. Similarly, we consider the averages of the product of three intensities given by 

G123 = (hhh) = (I 3 A ) + (4) + i(I 2 A )(h) + (II) (Ia)\ [3+2(cos^ 12 +cosfe+cos^ 3 )], (2.8) 
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and, assuming I a = h as before, we have 

(3) G$> a(3) a(2) / 3 \ 

Sl23 = ( h )(£}( h ) = V + 2 I 2 + cos ^2 + COS^3 + COS^i 3 J , (2.9) 

where g' 3 ' = (J^) / (Ia) 3 is the third-order normalized intensity autocorrelation function. 
In this paper, we consider the case of coherent light only. Then we have =g lyT> = 1. 

3 Event-by-event simulation 

We first discuss the general aspects of our event-by-event, particle-only simulation ap- 
proach. This approach is unconventional in that it does not require knowledge of the 
wave amplitudes obtained by first solving the wave mechanical problem nor do we first 
calculate the quantum potential (which requires the solution of the Schrodinger equa- 
tion) and then compute the Bohm trajectories of the particles. Instead, the detector clicks 
are generated event-by-event by locally causal, adaptive, classical dynamical systems. 
Our approach employs algorithms, that is we define processes, that contain a detailed 
specification of each individual event which cannot be derived from a wave theory such 
as QT. 

The simulation algorithms that we construct describe processes that are most eas- 
ily formulated in terms of events, messages, and units that process these events and 
messages. In a pictorial description, the photon is regarded as a messenger, carrying 
a message that represents its time-of-flight. In this pictorial description, we may speak 
of "photons" generating the detection events. However, these so-called photons, as we 
will call them in the following, are elements of a model or theory for the real laboratory 
experiment only. The only experimental facts are the settings of the various apparatuses 
and the detection events. 

The processing units mimic the role of the optical components in the experiment. A 
network of processing units represents the complete experimental setup. The standard 
processing units consist of an input stage, a transformation stage and an output stage. 
The input (output) stage may have several channels at (through) which messengers ar- 
rive (leave). Other processing units are simpler in the sense that the input stage is not 
necessary for the proper functioning of the device. A message is represented by a set 
of numbers, conventionally represented by a vector. As a messenger arrives at an in- 
put channel of a processing unit, the input stage updates its internal state, represented 
by a vector, and sends the message together with its internal state to the transformation 
stage that implements the operation of the particular device. Then, a new message is sent 
to the output stage which selects the output channel through which the messenger will 
leave the unit. At any given time, there is only one messenger being routed through the 
whole network. There is no direct communication between the messengers nor is there 
any communication between the processing units other than through the messengers. 
We view the simulation as a message-processing and message-passing process: It routes 
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messengers, representing the photons, through a network of message-processing units, 
representing the optical components in the laboratory experiment. From this general de- 
scription, it should already be clear that the process that is generated by the collective of 
classical dynamical systems is locally causal in Einstein's sense. 

3.1 Simulation model 

The network of processing units represents the whole experimental setup. For the present 
purposes, that is the demonstration that the HBT effect can be explained by a particle- 
only model, it is sufficient to simulate the bottom part of Fig. [TJ All the components, 
photons, beam splitters and photon detectors, have corresponding parts in our event- 
based simulation. As all the components are already presented in our previous work [32- 
35, 37-45], for completeness, we only give a brief description of each of the components 
of the simulation setup. 

3.1.1 Messenger 

We view each photon as a messenger. Each messenger has its own internal clock, the 
hand of which rotates with frequency /. When the messenger is created, the time of 
the clock is set to zero. As the messenger travels from one position in space to another, 
the clock encodes the time of flight t modulo the period 1 //. The message, the position 
of the clock's hand, is most conveniently represented by a two-dimensional unit vector 
e i = ( e o,/> e i,/) = (cosi/?/,sini/?/), where xpi =2nfti and the subscript I >0 labels the successive 
messages. The messenger travels with the speed of light c. In this paper, we do not need 
to specify the fixed frequency / and to specify a message, we use the angle ipi instead of 
the time-of-flight t\. 

3.1.2 Beam splitter 

The structure of the processing unit for a beam splitter (BS) is shown in Fig. |3l The unit 
has two input and two output channels labeled by k = 0,1 and consists of an input stage 
(DLM) a transformation stage (T), and an output stage (O). 

The input stage receives a message on either input channel or 1, never on both 
channels simultaneously. The input events are represented by the vectors v/ = (1,0) or 
v; = (0,1) if the Zth event occurred on channel or 1, respectively and are processed by 
a simple deterministic learning machine (DLM) [38-41,43]. The DLM has two internal 
registers Yjy = (Qy,Sjy) and one internal vector X; = (^o,/' x i,/)' where x ,/ + x i ,/ = 1 an d 
Xjy > for k = 0, 1 and all I > 0. Upon receiving the /th input event, the DLM performs the 
following steps: It stores the input message ejy = (cost/^ ^sint/^ j) in its internal register 
Y/ c ; = (Cjy,Sjy). Then, it updates its internal vector according to the rule 

x/ = 7*/-i + (l-7K (3-1) 
where < 7 < 1. A detailed analysis of the update rule Eq. can be found in Ref. [44]. 
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Figure 3: Diagram of a DLM-based processing unit that performs an event-based simulation of a beam splitter 
(BS). The processing unit consists of three stages: An input stage (DLM), a transformation stage (T) and an 
output stage (0). The solid lines represent the input and output channels of the BS. The dotted lines indicate 
the data flow within the BS. 



The transformation stage accepts the messages from the input stage, and transforms 
them into a new four-dimensional vector 



V2 



Q,/ a/^U ~~ So,/ ^/Xo} 



(3.2) 



The output stage sends out a messenger (representing a photon) carrying the message 

(3.3) 



w 



wo,/ 



where 



Wq,1 
si 



= (c ,i\Jxoj/2-S u ^x u /2^ I si, 
= ( C u yj x u /2 + S 0/l ^ x 0J /2j I s u 



(3.4) 
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Figure 4: Diagram of the event-based detector model defined by Eqs. (|3.7p and ()3.8p . The dotted line indicates 
the data flow within the processing unit. 



through output channel if sf > r where < r < 1 is a uniform pseudo-random number. 
Otherwise, if sf < r, the output stage sends through output channel 1 the message 



zo,i 
z l,i 



(3.5) 



where 



zq,i = [CijJ xy/2-Sm J Xqj/2) /ti, 



= [C 0r iJx 0i i/2+S u Jx u /2 ) ti, 



ti = \/ z l,i+ z li- (3-6) 

We use pseudo-random numbers to mimic the apparent unpredictability of the exper- 
imental data only: The use of pseudo-random numbers to select the output channel is 
not essential [39]. Note that in our simulation model there is no need to introduce the 
(quantum theoretical) concept of a vacuum field, a requirement in the quantum optical 
description of a BS. 

3.1.3 Photon detector 

A schematic diagram of the unit that functions as a single-photon detector is shown in 
Fig. H] [44]. The first stage consists of a DLM that receives on its input channel the Zth mes- 
sage represented by the two-dimensional vector e; = (cosipi,simpi) . In this paper, we use 
the simplest DLM containing a single two-dimensional internal vector with Euclidean 
norm less or equal than one. 

We write p; = {po,i,Pi,i) to denote the value of this vector after the /th message has 
been received. Upon receipt of the /th message the internal vector is updated according 
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to the rule 

P/ = 7P/-i + (l-7)e/, (3.7) 

where < 7 < 1 and I > 0. If 7 7^ 0, a machine that operates according to the update 
rule Eq. (|3.7|) has memory to store an amount of information that is equivalent to the 
information carried by a single mesasage only. Obviously, the rule Eq. \3.7) is the same 
as that used for the BS (see Eq. i3.1) ) but the input data is different. 

The second stage of the detector (see Fig. ID uses the information stored in the internal 
vector to decide whether or not to generate a click. As a highly simplified model for the 
bistable character of the real photodetector or photographic plate, we let the machine 
generate a binary output signal using the threshold function 

S, = ®(p?-r,), (3.8) 

where ©(.) is the unit step function and < T\ < 1 is a uniform pseudo-random number. 
Note that in contrast to experiment, in a simulation, we could register both the S/ = 
and S; = 1 events such that the number of input messages equals the sum of the S/ = 
and S; = 1 detection events. Since in experiment it cannot be known whether a photon 
has gone undetected, we discard the information about the S/ = detection events in our 
future analysis. The total detector count is defined as 

N = X>, (3.9) 

7=1 

where I is the number of messages received. Thus, N counts the number of one's gener- 
ated by the machine. 

3.1.4 Experiment 

The processing units that simulate the optical components are connected in such a way 
that the network corresponds to the experimental set up in the laboratory. As explained 
earlier, it is sufficient to consider the bottom part of Fig. [TJ 

4 Simulation results 

Following Ref. [59], the phase of the coherent photons emitted by the source is "random- 
ized" by letting the light pass through an EOM, the voltage of which is switched with a 
frequency of 50 Hz. To mimic this in the simulation, we send Ni nterm i messengers with 
some fixed but randomly chosen phase, then another Nj nterm i messengers with another 
fixed but randomly chosen phase, and so on. In practice, we use ^interval = 2500. The 
messengers (photons) are sent through either channel A or B, one at a time and are either 
transmitted or reflected by the beam splitters. Before hitting a detector, the messenger 
experiences a time delay corresponding to (pAn or (ps n (n = 1/2,3). The detector processes 
the message carried by the messenger and decides whether or not to produce a click. 
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Figure 5: Case 1: All BSs in Fig. [TJbottom) removed. Simulation results for the detector counts as a function 
of (p i = <pAi — tpgi . Top: The differences between the time-of-flights of the messengers entering channel A and 
the time-of-flights of the messengers entering channel B is constant. Bottom: The differences between the 
time-of-flights of the messengers entering channel A and the time-of-flights of th e m essengers entering channel 
B are random. Circles: simulation data; Dashed line: Wave theory solution Eq. (|2.2[l averaged over (p\. 



We consider three different experiments. In case 1, we remove all BSs in Fig.[lJbottom) 
and study the signal produced by detector D\ . Then, in case 2, we remove BS2, that is we 
consider the HBT experiment with two detectors, as indicated by the dashed-dotted line 
in Fig.[lJbottom). Finally, in case 3, we study the full three-photon correlation experiment, 
see Fig. [^bottom). I n cases 2 (3), the intensity-intensity correlations are calculated by 
counting coincidences of two (three) messengers, meaning that the arrival times of the 
two (three) messagers are within a time window W, to be discussed in Section 4.4. All 
simulations have been carried out with 7 = 0.99. 

4.1 Case 1: One detector 

Let us first demonstrate how the event-based model of the detector works [44]. The mes- 
sengers, randomly entering through channels A or B, are sent directly to the time-delay 
units that change the angle, representing the time-of-flight, by <pA\ or </>bi, respectively. 
The messengers are then processed by detector D\. We perform two different sets of 
simulations. First, we keep the differences between the time-of-flights of the messengers 
entering channel A and the time-of-flights of the messengers entering channel B con- 
stant. In this case, according to wave theory, we expect to see clear interference fringes. 
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Figure 6: Case 2: BS2 in Fig. [TJbottom) removed. Simulation results of the two-particle coincidence counts 
as a function of tp\i where (pyi = <p\ — (f>2, and <pn=<pAn~ < pBn ( w = l/2). The time-of-flights of the messengers 
entering channel A and the time-of-flights of the messen gers e ntering channel B are taken to be random. Circles: 
simulation data; Dashed line: wave theory solution Eq. (|2.7p . 



Second, the differences between the time-of-flights of the messengers entering channel 
A and the time-of-flights of the messengers entering channel B are taken to be random. 
Then, according to wave theory, there should be no sign of interference effects. As shown 
in Fig. |3 our particle-only simulation approach reproduces both features and the results 
are in very good agreement with the wave theoretical results (see Eq. (|2.2|) ). 



4.2 Case 2: Hanbury Brown-Twiss experiment 

We consider the HBT experiment with two detectors, that is we remove BS2 from the 
diagram in Fig. [^bottom). Messengers enter the apparatus through channel A or B, 
one by one. The time-of-flights of the messengers entering channel A and the time-of- 
flights of the messengers entering channel B are taken to be random hence, as shown 
in Fig. [^bottom) there is no first-order interference. When passing a BS, the message 
changes according to the rules explained in Section 3.2.1. Then, before entering the de- 
tector, the message is changed once more by cpAn or fgn ( n = 1/2), depending on which 
path the messenger took. If the two detectors fire with the time window W (see Sec- 
tion 4.4), we increase the number of coincidences. The simulation data shown in Fig. [6] 
confirm that this procedure reproduces the results of wave theory, see Eq. I|2.7|) . 



4.3 Case 3: Three-particle intensity-intensity correlation 

Finally, we consider the full correlation experiment Fig. QJbottom) with three detectors. 
The simulation procedure is the same as in case 2, except that we count coincidences of 
clicks of three different detectors. Also in this case, the simulation data shown in Fig. [7| 
confirm that this procedure reproduces the results of wave theory, see Eq. I|2.9|) . 
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Figure 7: Case 3: Three particle correlation experiment (see Fig. [TJbottom)). Simulation results of the three- 
particle coincidence counts as a function of tp\i where <p\2 = <p\ — <p2, ar| d <pn=<pAn~ ( pBn ( w = l>2,3). We only 
show data for the case <pA2 = ( pB2 = Q' <Pa\ =( Pb?>< <Pbi =( Pa3 where (p^i and (p^\ are chosen randomly. The 
time-of-flights of the messengers entering channel A and the time-of-flights of the messenge rs en tering channel 
B are taken to be random. Circles: simulation data; Dashed line: wave theory solution Eq. I|2.9p . 



4.4 Discussion 

Our simulation model is based on a particle picture and makes no reference to concepts 
or results from wave theory. In contrast to the conventional quantum theoretical expla- 
nation in terms of the wave-particle nature of photons, our simulation approach requires 
a particle picture of photons only. During the event-by-event simulation we always have 
full which-way information of the photons (messengers) since we can always track them. 
Nevertheless, depending on the settings of the optical apparatuses, intensity-intensity in- 
terference is observed. Although the appearance of an interference pattern is commonly 
considered to be characteristic of a wave, we have demonstrated that, as in experiment, 
it can also appear as a result of a collection of particles that interact with the various op- 
tically active devices such as beam splitters and detectors. In this paper, we considered 
the case that is equivalent to a light source that produces photons in a coherent state only. 
The case of a thermal light source will be considered in future work. 

In real experiments, and also in our simulation approach, it is necessary to specify the 
procedure by which we count coincidences of detection events. For the experiments at 
hand, one introduces a time window W and one defines as a two (three) particle coin- 
cidence, two (three) detection events with the time difference(s) are smaller than W. As 
discussed extensively in our work on the simulation of Einstein-Podolsky-Rosen-Bohm 
(EPRB) experiments [35], the choice of the time window W is of crucial importance, both 
in the simulation and in real experiments [48], to obtain the correlation of a quantum 
system in the singlet state. In general, only when W — > 0, experiment and simulation can 
reproduce the correlation of a quantum system in the singlet state [35]. For large enough 
W, the relation to a quantum system in the singlet state is lost. In this paper, we have cho- 
sen W sufficiently large and generated groups of two (three) messengers such that if the 
two (three) detectors fire, this constitutes a coincidence of two (three) particles. In other 
words, the time delays are only used by the detector but are ignored in determining coin- 
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cidences. In this sense, the simulation results presented in this paper pertain to classical 
light and are therefore in excellent agreement with classical wave theory. To study the 
quantum aspects of two- and three-particle correlations the time delays should be used 
to also determine the coincidences, as in our EPRB simulations [35]. We leave this very 
interesting topic for future research. 

5 Conclusion 

We have demonstrated that our classical, locally causal, particle-like simulation approach 
reproduces the results of the Hanbury Brown-Twiss effect. Our event-based simulation 
model, a classical, locally causal, adaptive dynamical system, reproduces the results of 
wave theory without making reference to the solution of a wave equation and provides 
a simple, particle-based mental picture for what each individual photon experiences as it 
travels from the source to the detector. Our simulation algorithm demonstrates that the 
wave-particle duality is not the only way to describe the nature of a phtone but that there 
is another way that only needs the particle nature, satisfies Einstein's local causality and 
does not defy the common sense. Finally, we would like to emphasize that the algorithms 
used to simulate the optical components in this paper have not been designed to simulate 
the HBT-type experiments. The algorithms have been taken, without modification, from 
our earlier work on very different quantum optics experiments [32-35,37-44,44,45]. In 
this sense, it seems that our approach has predictive power: The algorithms can be reused 
to simulate very different experiments than those for which they were originally devel- 
oped. 
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